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.  INTRODUCTION 


In  recent  years  the  need  has  arisen  at  ARL  from  time  to  time  for  the  calculation  of  rigid- 
body  gust  response  characteristics  of  aircraft  configurations.  Such  calculations  are  useful,  for 
example,  in  assessing  the  platform  stability  of  air  vehicles  in  a  turbulent  atmosphere.  Apphcali.  i.s 
include  comparisons  between  competing  RPV  designs  and  design  requirements  for  a  t.  reel 
designator.  Similarly  ,  estimates  can  be  made  for  vehicle  flight  path  dispersion  or  for  gust  k  ads 
on  different  parts  of  the  vehicle.  Further,  Military  Specifications  (e.g.  Ref.  I )  call  for  calculation 
of  aircraft  response  to  atmospheric  disturbances  for  assessment  of  flying  qualities. 

This  Note  documents  a  general  gust  response  program  which  has  been  developed  at  ARL 
in  order  to  meet  such  requirements.  Although  the  program  has  been  configured  to  produce 
response  spectra  due  to  isotropic  homogeneous  turbulence  with  a  von  Karman  spectrum  it  is 
readily  adaptable  to  more  complex  turbulence  models  including  cross-correlations.  On  the  other 
hand  the  program  may  simply  be  used  to  provide  aircraft  gust  response  transfer  functions  which 
may  he  used  in  calculating  the  response  to  a  discrete  gust.  Most  of  the  parameters  describing 
the  aircraft  configuration,  llight  conditions  and  aerodynamics  are  read  in  as  inputs  thus  facili¬ 
tating  comparisons  between  different  aircraft  configurations,  investigation  of  effects  of  changes 
in  aerodynamic  parameters,  control  system  operation,  etc.  Although  limitations  exist  in  the  range 
of  validity  of  the  theory  it  may  be  assumed  that  its  use  in  assessing  relative  trends  extends  over  a 
rather  wider  range. 

Section  2  of  this  Note  summarizes  the  background  theory  and  treats  in  more  detail  the  air¬ 
craft  and  gust  mathematical  models  including  a  discussion  of  the  range  of  validity  implied  by 
the  various  approximations.  This  is  followed  in  Section  3  by  a  description  of  the  computer  pro¬ 
gram  and  Section  4  provides  examples  of  its  application  to  two  widely  different  vehicle 
configurations. 


2.  THEORY 

The  theoretical  treatment  in  this  section  is  derived  largely  from  References  2  and  4.  The  basic 
mathematical  approach  is  first  described  followed  by  a  more  detailed  look  at  the  two  basic 
components,  viz.  the  aircraft  dynamics  model  and  the  treatment  of  the  gust  input. 

2.1  Mathematical  Model 

The  aircraft  response  to  gust  input  can  be  treated  as  a  small  disturbance  problem  about  a 
nominal  trim  state.  Thus  a  linearised  set  of  equations  provides  an  adequate  description  of  the 
motion.  A  general  Laplace  Transformed  set  of  equations  can  be  written  as  follows: 

Afv)x(.v)  =  f(x)  =  B(.v)g(.v).  (I) 

Where  x(.s)  is  the  Laplace  Transform  of  the  aircraft  state  vector.  f(.v)  is  the  Laplace  Trans¬ 
form  ol  the  aerodynamic  force  vector  due  to  the  gusts  and  the  matrix  A(.v)  relating  the  two  con¬ 
tains  a  description  of  the  aircraft  aerodynamic  and  inertial  characteristics.  The  force  vector. 
f(.v),  can  be  derived  from  a  representation  of  the  gust  velocity  field  g(.v)  through  B(.v).  which  is 
a  matrix  of  "gust  transfer  functions".  These  w  ill  be  treated  in  more  detail  in  2.3  below  . 

From  Equation  (I)  a  matrix  of  Transfer  Functions,  G(.v).  relating  the  aircraft  response. 
x(.r),  to  the  gust  input,  g(.v),  can  he  obtained: 

x(.v)  =  G(.v)g(.v).  (2) 

If  x(.v)  is  an  (//  x  I)  vector  and  g(.v)  is  an  (m  x  I)  vector  it  follows  that  G(.v)  is  an  (»  x  m) 
matrix.  A  substitution  of  iw  for  s  (<«  being  frequency  in  radians/s)  yields  the  Fourier  Transform 


representation  of  the  transfer  function  G(iw).  If  the  response  to  a  discrete  gust  is  required  then 
this  can  be  done  directly  in  the  frequency  domain  through  the  use  of  G (iw)  followed  by  a  trans¬ 
formation  back  to  the  time  domain.  On  the  other  hand,  a  spectral  representation  of  the  response 
can  be  obtained  immediately  from  the  relation: 

4v  =  C**gGr.  (3) 

Where  the  matrix  is  gust  input  power  spectral  density  matrix  tin  x  in)  and  4v  is  the  out¬ 
put  power  spectral  density  matrix  (n  x  /;).  Normally  the  diagonal  elements  of  4V  are  of  interest 
and.  for  these.  Equation  (3)  can  be  expanded  as  follows: 

AV.-t.I***)  =51  X  G*k{ia>)'frllkfi1(w)G,l(iw)  (4) 

*i/i 

where  the  subscripts  refer  to  elements  of  the  respective  matrices.  For  the  case  of  zero  input 
cross  spectra,  i.c.  no  cross-correlations.  Equation  (4)  can  be  further  simplified  to  give 

m 

'IV,.r,(ud  =  X  >)  to).  (S) 

*  1 

The  integral  of  the  output  power  spectrum  with  respect  to  w  gives  the  variance  of  the  air¬ 
craft  response  about  its  nominal  steady  value,  i.c. 

a  V,  =  |*<IVIjJ(<t»)(/u/.  (5a) 

where  <V,  is  the  variance  of  the  .v,  component  of  the  state  vector. 


2.2  Aircraft  Model 

The  use  of  linearised  small  disturbance  equations  allows  the  aircraft  dynamics  to  be 
decoupled  into  Longitudinal  and  Lateral  response  sets,  thereby  reducing  the  size  of  the  A(.v) 
matrix  (Eqn  I).  The  equations  can  be  augmented  with  additional  equations  describing  any 
flight  control  systems  that  may  be  present. 


2.2.1  Aircraft  Dynamics 

The  decoupled,  non-dimensional,  small-disturbance  equations  used  arc  taken  from 
Reference  2.  The  nominal  reference  state  is.  in  general,  steady  flight  at  a  climb  angle  of  y0.  Tile- 
thrust  vector  is  assumed  to  be  at  an  angle  w  ith  respect  to  a  reference  set  of  axes  fixed  in  the 
body.  This  set  of  reference  axes  is  assumed  to  be  initially  aligned  with  the  velocity  vector,  i.e. 
stability  axes. 

The  Longitudinal  equations  arc  in  wind  axes  with  the  state  vector  given  by: 

x(.r)  =  [AF(.s),  ot(.s),  q(s).  Ad(x)]7'  (6) 

while  the  lateral  equations  are  in  body  axes  and  the  state  vector  is: 

x(.c)  =  [/3(s).  p(s).  /-(*-).  rf>(s  )] r.  (7) 

The  A(.r)  matrix  for  each  of  these  cases  is  given  in  Appendix  I.  Note  that  since  non- 
dimensional  equations  are  used,  the  Laplace  Transform  is  with  respect  to  non-dimensional  time. 
(2F/c)r.  Further,  for  generality,  the  aerodynamics  characteristics  are  written  as  aerodynamic 
transfer  functions  even  though  in  practice  these  are  replaced  by  their  quasi-steady  approximations 
in  terms  of  aerodynamic  derivatives. 

Although  Equations  (6)  and  (7)  list  only  the  state  variables,  responses  for  related  variables 
can  readily  be  derived.  For  example  the  variation  in  load  factor.  Sn(s)  in  g's.  in  response  to  a 
vertical  gust.  u„  say.  is  related  to  a(.v)  and  the  incidence  and  airspeed  responses  due  to  the 

vertical  gust,  respectively: 

Ms)  =  i.*(r))  I  2 A  F(x).  (8) 


2 


Similarly,  the  yaw  angle  Ms)  is  obtained  from  the  integral  of  the  yaw  rate  /(w  !>y  : 

{ .  R )  .  ,v  .  M  v )  see  y,  .  r(  s)  ( 9 1 

where  -R  is  the  ratio  of  span  to  mean  aerodynamic  chord. 


2.2.2  C  ontrol  Systems 

The  etfeet  of  automatic  llight  control  systems  on  the  gust  response  can  he  obtained  hy 
extending  the  state  vectors,  liquations  (6)  and  (7).  and  the  At s i  matrix  to  model  the  systems 
present,  for  example  the  pitch  or  yaw  damper  systems  on  the  Mirage  III  may  he  modelled  hy 
the  following  equation: 

t  '  i  v  ■  I )( r.;v  ■  litra.t  I  )-■■>(  V>  -  A.s.i/lvl  t  l<  >) 

where  <>t .v >  is  the  damper  response  to  pitch  rate,  i/i'i.  A.  is  a  constant  gain  (depending  on  speed 
and  height)  and  n.  rj  are  time  constants,  ltnis.  for  example,  the  longitudinal  state  vector, 
fquation  < r> ).  need  only  he  augmented  hy  >(v)  while  the  A(v)  matrix  has  an  additional  row. 
[0.0.  Av.  0.  /hv)].  and  an  additional  column.  {  tV i .  -  (m  .  ( /,„  .  0.  Ii(s)\r.  where  /i(v)  - 
( "|.v  litres  l)(r-.t.v  ■  I  I-. 

In  practice,  the  computer  program  described  in  Section  a  allows  only  quadratic  functions 
of  s  in  the  A(vi  and  Bts)  matrices.  I  his  is  readily  accommodated  hy  introducing  an  additional 
state  variable,  mf')  sav  and  replacing  fquation  1 10)  with  the  two  equations: 

(r.sV  1  1(7.0  |  |.S(V|  -  n„A)  (  (  ( 

(r.s  I  )l  r.o  l|.',.l»l  Ax/ |v) 

The  equivalent  Aid  matrix  is  shown  in  .Appendix  I. 


2.3  Ci list  Model 

I  his  section  deals  with  the  right-hand  side  of  fquation  I.  Different  levels  of  approximation 
for  the  gti't  input  are  first  described  followed  hy  a  discussion  of  the  gust  spectra  (o  he  used.  It  is 
assumed  throughout  that  the  aircraft  samples  a  frozen  gust  field  on  a  straight  lire.  i.e.  there  is  no 
departure  of  the  aircraft  from  rectilinear  flight . 


2.3.1  Point  Approximation 

1  he  gust  field  is  characterized  hy  the  three  turbulent  velocity  components  i/L..  ig.  and  irc.  In 
the  point  approximation  the  variations  of  these  velocities  over  the  aircraft  are  ignored,  i.e.  the 
aircraft  is  treated  effectively  as  a  point.  Consistent  with  the  small  disturbance  assumptions  the 
gust  vector,  g.  is  separated  into  longitudinal  and  lateral  sets.  j>i  and  gj  respectively.  Non- 
dimen'ionahsing  hy  dividing  hy  I',,  and  noting  that  in  the  point  approximation  ug  I,. 
and  r,;  I  ,.  -  tiK.  the  two  sets  become: 


gi  I  ,-.  H-Ir 

Si  -  If’cl r 


It  follows  that  the  B(.v)  matrices  (fqn  1 )  required  to  relate  these  to  aerodynamic  forces  are: 


B,(.v)  = 


f'liv  (fTV  (/lij, 

f'l.v  (/|., 

(>  i,i 

:  Hats)  - 

hmV  Cm  T 

('n/; 

0  (1 

0 

(13) 


The  aerodynamic  transfer  functions,  (i'tv.  (/nv.  etc.,  are  usually  replaced  with  their  quasi- 
steady  approximations  (see  Appendix  If 


3 


2.3.2  Power  Series  Method 


Phis  is  an  extension  to  the  point  approximation  in  that  the  east  t  eld  is  now  represented  by 
the  gust  xelocities  and  their  first  derixatixes.  It  is  still  assumed  that  variations  m  the  c-dircction 
(normal  to  the  flight  direction)  are  negligible.  1  he  remaining  two-dimensional  velocity  held  is 
expanded  as  a  Tax  lor  series  about  the  aircraft  centre  ot  gravity.  and  higher  order  terms  are 
neglected,  e.g. 

«>K(  v.  t  . /)  -  «»,.(/)  'fix.  x  mi,,  i  t  .1  higher  order  terms.  I  I4i 

l  or  aircraft  rigid-body  responses  Reference  1  suggests  that  it  is  suflicient  to  consider  only 
the  uniform  gust  immersion  terms  m.  rB.  nf  together  xxith  linear  gradient  terms 

(I  I  ,1  .  <  II.-  I  X  l/j 

1 1  I  ,.)  .  <  n  •  .  i  />„.  t  i 

(I  I,  )  .  1 1  c  <  ,\  rc 

Further.  i/t..  (■„.  n  ^  and  can  be  considered  to  be  mutually  independent  ( uncorrelalcd )  xx  hile 
correlations  haxe  to  be  accounted  for  between  n.  and  r/k-.  and  between  rc  and  r(.  Now  Reference  3 
notes  that,  of  the  gradient  terms,  /c  is  dominant;  consequently  onlx  this  term  is  retained  in  the 
folloxxing  treatment,  for  simplicity. 

The  longitudinal  gust  set  thus  remains  the  same  as  in  I  cjuution  (12)  while  the  lateral  set 
becomes : 

S.  ==  I  4-.  p,-l r  1 1  hi 

where  the  /c.  gust  can  be  thought  of  as  corresponding  to  a  roll  rate  I  he  cx|iii\aleni  B.ni  matrix 
(bun  (13))  becomes: 


2.3.3  Range  of  Validity 

Assuming  rectilinear  flight  at  speed  I  ,,  and  neglecting  variations  in  the  ^-direction  the  two- 
dimensional  gust  spectral  component  as  experienced  by  the  aircraft  is  a  xelocitx  held  of  the  form 

(.„U| U..1/1  (|S| 

which  represents  a  time  periodic  xelocitx  at  any  point  !  v.  r)  of  the  \ chicle.  Ui  and  LJj  are  related 
to  the  waxelengths  A,  and  A.-.  in  the  x  and  i  directions,  i.e. 

o,  2tt  (Il>) 

The  range  of  validity  of  the  point  and  power  series  < or  linear  held)  approximations  can  be 
stated  in  terms  of  1>,  and  IT.  Reference  5  notes  that  the  linear  neld  approximation  goes  a  good 
representation  of  the  gust  velocity  distribution  oxer  the  aircraft  only  when  i>,  and  are  “small 
enough",  i.e.  when  A,  and  A_.  are  considerably  (e.g.  ten  times)  larger  than  the  length  and  span, 
respectively .  of  the  aircraft. 

Further,  considerations  from  unsteady  aerodynamics  provide  limits  on  11,  and  U  -  for  which 
the  quasi-steady  approximation  is  acceptable.  Thus  the  following  approximate  limits  arc- 
obtained  (Refs  2,  5): 

LI,  .  ),<•  <  0  1, 
i> . .  />  ^  2. 
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(  20) 
(21) 
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These  are  roughly  compatible  with  the  limitations  on  t|.  a  due  to  the  linear  held  represen¬ 
tation.  implying  a  consistent  use  of  quasi-steady  acrody naniies  with  i lie  linear  held  model.  Ir. 
general,  the  range  ot  frequencies  covered  by  such  a  treatment  is  adequate  tor  rigid  hods  mainin-. 
appropriate  to  stability,  guidance  and  control  problems,  but  appreciable  errois  ean  isoir  d 
structural  mode  responses  are  important,  l  or  this  case  a  more  aecurate  act.idyiumic  theory 
should  be  employed  and.  consistent  with  this,  a  more  exact  sp"  al  disiribution  ol  the  gu-t  held 
(see.  for  example.  Ref.  5). 


2.3.4  Gust  Spectra 

I  he  simplest  model  of  turbulence  assumes  that  above  the  planetar,  boundary  laver 
(approximated  2lHK)ft)  the  turbulence  is  isotropic  and  homogeneous.  I  his  means  that  length 
scales  and  intensities  are  the  same  in  all  directions.  I  he  von  Kantian  spectra  lor  the  turbulence 
velocities  are  gi veil  by  (Ref.  I)' 

,h„,.„,(,-!i)  ",-(;/  -i  ;  I  tl  'W/o, t;:i 

«l>  ,(ii|  I  il-hi  /  r  i .  ;  I  -  a  I  ■  .v'V/ LJ|  i-;  I  1 1  •  339/ lit)-:1’  \  (23l 

where  \  ■  i-  ti  e  turbulence  intensity  and  /  the  length  sc.de.  taken  to  be  25otMt  Itoni  Relereiive  i. 
■Ml  ve'ocity  cross  spectra  are  zero  in  isotropa.  turbulence.  However,  cross  spectia  arise  even  m 
.sotropic  turbulence  if  variatii'iis  in  velocities  over  the  vehicle  arc  allowed  isec.  h'r  example. 
Kef  -I/.  ( liese  are  neglected  here  as  explained  in  Section  2. 2.  I  he  only  other  spectrum  required 
here  is  ttiat  for  given  by  Reference  I  • 


<1>,  .^.tU, i  -  t-v  /  tit  n -/  -1  hi- 


1 24 1 


where  t<  is  the  span. 

All  spectra  above  are  single  sided  so  that  the  integration  limits  in  I  qualionf  5a) are  from  zero 
to  mluuty .  Note  that  il.  and  frequency.  <•>.  arc  related  by  i.\  I  ,  (see  I  q n  i  IX)). 

Mthough  a  simple  isotropic  model  of  turbulence  has  been  described,  the  general  approach 
used  makes  it  relatively  simple  to  use  a  more  complex  niodci  appropriate  tv',  say.  (light  at  lower 
altitudes. 


3.  COM  PI  "I  HR  PROGRAM 

I  he  mam  program  has  been  designed  to  be  as  flexible  as  possible  to  enable  the  user  to  alter 
it  to  suit  the  particular  problem  in  hand.  I  luis  the  mam  program  concentrates  on  defining  the 
problem  parameters  while  the  basic  calculations  arc  achieved  by  calls  to  general  purpose  sub¬ 
routines.  It  is  also  possible  to  use  the  program  simply  to  generate  transfer  functions  for  any 
system  ol  the  form  given  by  liquation  I .  A  general  description  of  the  program  follow  s  including  a 
summary  of  the  input  requirements  and  output  produced  by  the  current  gust  response  calcula¬ 
tions.  I  lie  program  is  written  in  I-ortran  and  implemented  on  the  ARL  I) PC'  System  10. 

3.1  General  Description 

The  basic  structure  of  the  main  program  GISTR  is  shown  in  figure  I.  The  problem 
dimensions  ti.e.  dimension  of  the  A  matrix)  and  all  variable  dimensions  used  in  (he  program 
are  defined  at  the  beginning  of  the  program  in  a  Parameter  statement.  This  is  followed  by  a 
block  with  labelling  and  initialisation  data  required  for  later  plotting  of  the  results.  The  plot 
tiles  for  the  Calcomp  plotter  arc  prepared  via  the  program  “TRANS"  described  in  Reference  6. 

The  next  step,  constituting  the  problem  definition,  is  to  specify  the  elements  of  the  Atxi 
and  B(.v)  matrices  (I  qn  ( I )).  Each  element  of  A(.v)  and  B(.x)  is  assumed  to  be  a  second  order  poly¬ 
nomial  in  v.  i.c. 

(h  i  -  «i''  i  flav¬ 
or  (25) 

h»  1  h\s  :  /C..V-. 
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In  order  to  define  Af.v).  three  arrays  ol  coctlicients  .ire  set  up.  the  first  to  eont.nn  all  the  zero 
order  coeflicients  («,*.  ete.).  the  second  to  contain  the  first  order  coefficients  t,;,.  eic.  i  and  the  third 
to  contain  the  second  order  coefficients  («■_>.  etc. I  Ihree  tuo-diniertsion.il  aruvs.  -fO( l.Jt. 
dll/../)  and  (2(  /.  ./  )  are  used  for  this  purpose  such  that  lOt  I.  ./  i  contains  the  a,,  coclhcient  ot 
the  (/../)  element  of  A(v)  etc.  Similar  considerations  apple  to  lini  except  (fiat,  since  transfer 
functions  are  calculated  for  one  gust  component  at  a  time,  one  dimensional  arrays  H0[!  H  l(/ 1 
and  />’2l / )  are  required. 

I  or  the  aircraft  case  the  elements  c>f  the  Al'l  matrices  Ii.iu-  been  listed  m  Appendix  1  ,.nd 
the  equivalent  Itfr  I  matrices  arc  given  by  I  c|ti.itioiis  1 1  o  and  ( I  '). 

Having  defined  the  system  matrices  the  solution  proceeds  with  a  call  to  subroutine 
"7  R  ANSI  ".  which  calculates  the  polvnomi.il  coellicients  ol  the  denominat«>r  and  iH.mer.itois 
of  each  of  the  transfer  functions  defined  by  the  svstem.  I  tie  detiomm.itoi.  or  characteristic. 
polvnomi.il.  root'  are  obtained  hv  a  call  to  subroutine  "1*01  R  I  .  I'olv  iiomi.d  coethcient s 
of  any  extra  transfei  functions  requited  can  be  calculated  at  this  stage.  I  he  present  program 
calculates  those  for  the  normal  acceleration  at  the  centre  of  gravity  t!  qnisn  in  the  I. oneiiiidm.il 
case  and  the  yaw  ancle  (I  qn  iRt)  in  the  l  ateral  case. 

I  he  treipienex  response  tunctions.  (i  ;.u.. o.  are  obtained  by  substituting  v  n  ■:  2. 

where  the  rion-dimensionahsing  time,  r  f2l  ,.  appears  because  the  l.apiace  variable  v  is 
.Mth  respect  to  non-dimensional  time.  I  he  response  power  spedial  densities  then  follow 
from  f  quabon  Oi  using  the  spectra  eiven  in  Section  2.3  4.  In  making  these  calculations  use  :s 
m  ide  of  sii hi i>u line  "C'I’Ol."  w  Inch  ev  aluates  a  poly  noun. 1 1  h.iunti  complex  coeliic icnts.  '(  l’(  )L" 
is  appended  to  the  main  procram. 

Computation  nine  is  determined  largely  hv  subroutine  I  R  ANSI".  If  V»i  is  <i  ■  c  then 
"  !  R  \NSI  '  has  to  evaluate  !  >  determinants,  f  I. to  CIM  time  is  ptopertion.il  to  •:  1  i. 

l  or  /:  i'  t ins  amounts  t c*  approximately  (t  seconds  on  the  1)1  C  System  10.  increasing  verv 
rapidly  for  larger  If  larger  order  systems  need  to  be  treated,  computation  time  can  he  reduced 
hv  allow  ini:  ft  pc  her  order  polynomials  ;i'  elements  of  A/.vt  and  Jhs).  fhis  would  require  an 
alteration  to  subroutine  "  I  K  ANSI  ". 

3.2  Inputs 

l  lte  data  required  to  define  the  matrices  At  t  >  and  Hi  \  i  w  ill  depend  on  the  particular  problem 
m  hand.  In  general  the  lata  can  he  divided  into  three  croups  related  to  aircraft  conlicuialion. 
Ihclit  conditions  and  aerodynamics.  As  an  example  Appendix  2  lists  the  set  of  data  required  to 
calculate  longitudinal  and  lateral  gust  responses  of  the  Mirage.  I  he  pitching  and  yawing  moment 
data  are  oltcn  given  relative  to  a  reference  centie  of  gravity  position  and  the  program  adjusts 
these  to  the  actual  c.c.  t  he  drag  due  to  pitch  damper  derivative.  G>  .  has  been  sei  to  zero  to  avoid 
a  non-linearity  arising  from  its  sign  depending  on  lire  sign  ot  a  If  is  in  any  ease  small.  It  has  also 
been  assumed  that  the  reference  state  is  one  of  level  flight.  :  ,  0.  and  that  the  thrust  line  is 

parallel  to  the  reference  axes.  i.e.  r-|  -  0. 

With  the  given  data,  the  program  carries  out  all  the  necessarv  non-dimeti'ionahsing  and 
auxth.u  v  calculations  necessarv  to  set  the  matrix  elements.  One  additional  piece  of  information 
required  is  the  thrust  due  to  speed  derivative,  G,  f  or  a  jet  powered  aircraft  such  .is  the  Mirage, 
constant  thrust  propulsion  can  be  assumed  so  that  Or,.  -  2 On 

Having  set  the  elements  of  A(.r)  and  B( v )  and  calculated  the  transfer  function  elements. 
(Via.  the  program  only  needs  the  turbulence  length  scale  to  define  the  input  spectra  ne-  (i.e.  from 
I  i|ns  (22)  (24))  and  obtain  the  corresponding  output  spectra  <r,.-  from  l.quation  (3).  where  a,.- 
is  the  relevant  gust  variance. 

3.3  Ouputs 

The  choice  of  output  depends  on  the  application  and  can  be  readily  changed  by  the  user. 
Immediately  available  are  the  coeflicients  of  llte  Laplace  form  of  the  transfer  functions.  Cl, 
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the  characteristic  roots  of  the  system.  the  real  anti  imaginary  parts  of  the  transfer  functions  as 
functions  of  frequency.  </,*(! *i ).  anil  output  spectra  as  Itinclmn s  i»f  fretpienc).  i-h  <«  I 

l  or  plots  of  frequency  response  or  output  spectrum  against  frequency  it  is  getter ally  con¬ 
venient  to  use  a  log  scale  for  the  abscissa  because  of  the  wide  liequency  range  generally  of 
interest.  In  this  case  it  is  useful  to  plot  spectral  density  multiplied  In  frequency.  i.e.  Upl'ilii  I 
rather  than  <I>(U|).  as  the  ordinate  since  the  area  under  such  a  graph  is  then  proport-onal  to  the 
variance  (1  qn  (5a)): 

I  iiplHUi  )(/(logmUi )  -  I  Upl'd  »i)</(  In  U,)  In  It)  (I  In  lull  >l,(12l  )</<2i  .1-' In  10 

!  0  J  0  J  II 

rims 

2  5  I  lipl'f. 'i  li/i loguiUi  I.  (2b) 

I  o 


4.  I  WMIM  I  S 

I'vpicai  results  arc  shown  lor  two  widely  dillcring  configurations  in  this  section.  I  he  lii'st 
one  can  be  thought  of  as  ;i  possible  RP\  configuration  while  the  second  one  relate"  to  the 
Mirage  III.  I  he  usefulness  of  the  program  in  assessing  the  effects  on  gust  response  of  changes  in 
c  g.  position,  control  'Vstems  and  aerodynamic  parameters  is  demonstrated. 

I  wo  graphs  are  plotted  for  each  response  variable  of  interest.  I  he  upper  graph  plots  login 
of  the  square  of  the  Dequeues  response  and  the  lower  graph  plots  (heoutput  power  spectrum  in  the 
form  llpbU-h  i  w  here  •»*-'  is  the  appropriate  gust  variance.  In  each  case  the  abscissa  is  logiutD i ) 
where  I  .  1 1  m  ft  s  here).  Since  the  present  program  defines  the  input  spectra  to  be  half 

those  given  by  fquations  (22)  (24)  the  output  variances  are  4  6  times  the  area  under  the  output 
spectra  graphs  (I  qn  (2f>)). 


4.1  RI’V  Configuration 

Results  shown  in  figures  2  4  are  for  a  vehicle  of  chord  0  I  I  m  (0  56  ft)  and  aspect  ratio 
10  living  at  106  m  s  ( 54S  ft  s)  at  an  altitude  of  505  m  ( 10(H)  In.  for  this  case  different  turbulence 
length  scales  are  used  for  the  vertical  gusts  (/  505  m  in  uu-  and  pK  spectra)  and  the  longitudinal 

and  lateral  gusts  (/  561  m  in  ty  and  ty  spectra). 

f  igure  2  shows  the  pitch  rate.  </.  response  to  vertical  gust.  «■„.  as  a  function  of  static  margin. 
As  the  c.g.  moves  forward  the  static  margin  increases  from  a  low  value  of  0  I  to  a  maximum  of 
I  -47.  At  the  same  lime  the  vehicle  becomes  much  more  sensitive  to  vertical  gusts  Ihe  range  of 
validity  of  the  calculations  (I  qn  (20))  is  also  indicated  on  figure  2.  I  qinvalent  results  for  the 
pitch  attitude.  0,  response  are  shown  in  figure  5.  Ihe  two  obvious  sets  of  peaks  correspond  to 
the  short  period  and  Phugoid  frequencies  of  the  vehicle,  f  or  static  margin  of  0  I  there  is  very 
little  pitch  attitude  response  for  logml'i  above  2  (i.e.  «<  above  5-4S  rad  s)  in  contrast  to  a  con¬ 
siderable  amount  of  energv  between  2  log  U[  v.  I  5  (5-4S  »-  <,i  •.  II  (X)  rad  s)  for  a  static 
margin  of  I  -47.  In  this  latter  case  the  area  under  the  curve  between  these  limits  is.  verv  approxi¬ 
mately.  10  (square  inches)  x  5  x  10  7  (scale)  which  contributes  to  a  variance  «- •nt- 
4  6  x  5  x  10  7  2  5  x  10  B  (rad -/(ft  s  ')-).  for  a  gust  intensity.  <w-  of  10  ft  s  this  would 

mean  a  pitch  attitude  response  standard  deviation  of  0  01 5  rad  0  X7  degrees  in  the  given 
frequency  range. 

An  example  of  roll  rate.  />.  response  to  lateral  gust,  rK,  and  gust  gradient.  pK.  is  shown 
in  figure  4.  The  peaks  near  log(12|)  2  correspond  to  the  Dutch  Roll  frequency.  The  roll 

rate  gust  response  variances  due  to  ly.  and  /y  can  be  obtained  from  the  areas  under  the  respective 
output  spectra  over  the  desired  frequency  range,  as  was  done  above  with  the  pilch  attitude  case. 
Since  it  has  been  assumed  that  there  is  no  cross  correlation  between  the  t\,  ami  /v  gusts  the  total 
output  variance  for  roll  rate  is  given  by  the  sum  of  the  variances  due  to  and  to  /y 
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4.2  Mirage-  III  Configuration 

I  tgtircs  5  S  'how  results  lor  roll  rale  and  yaw  rate  responses  ol  a  Mirage  tspe  configuration 
to  lateral  gusts.  r£.  These  are  lor  flight  at  M  =  0-9  at  6lt)  m  ( 2tKM)  ft)  altitude.  I  he  turbulence 
length  scale  has  been  taken  as  2500  ft. 

figures  5  and  (>  show  the  influence  of  dihedral  e fleet,  (T  .  on  the  roll  rate  and  yaw  rate- 
responses  respectively.  The  Dutch  Roll  frequence  (about  4  rail  s).  and  the  roll  mode  time¬ 
s'  instant  (about  0  25  sido  not  depend  strongly  on  (i..  but  the  spiral  mode  time  constant  changes 
from  about  400  s  for  Ci  0-l)l5(closc  to  the  true  value)  to  a  sets  low  (is  foi  Ct  0  15. 

I  lie  effect  on  the  roll  response  to  lateral  gust  is  also  scry  large,  as  mat  be  expected  (I  ig.  5).  but 
the  saw  rate  response  only  changes  slightly. 

I  null',  the  effect  of  saw  damper  on  roll  and  saw  rate  responses  is  shown  in  I  igures  ?  and  V 
Ihe  damper  is  clearly  effect  is  e  around  the  Dutch  Roll  frequencs  and  consequentls  is  also 
effect  is  c  m  diminishing  the  roll  and  saw  gust  responses  ss  liiclt  are  pi  edom  mantis  at  this  frequencs . 
Ihe  approximate  range  of  validity  shown  in  I  igure  X  covets  most  ol  the  output  spectium  esen 
for  tills  case  of  large  chord.  I  lie  range  of  validity  on  li i  can  be  extended  by  a  (actor  of  10  if 
unsteady  oscillatory  aerodynamics  are  used. 

Similar  plots  can  readily  be  piiHiuced  lor  ,m>  other  sariables  oi  interest  including  alls  of 
ills'  'la'.'  variables  (I  s|iis  i(i)  and  ('ll  or  extra  sariables  sus.lt  a-  itonnal  acceleration,  sass  angle, 
etc.  te.g.  I  s|iis  iS)  and  (9)i.  II  mils  the  responses  to  the  gust  components  gisen  bs  I  quations 
i!2t  ami  (i(0  are  considered  then  cross-correlations  can  be  neglected  and  I  quation  (5l  can  be 
used  to  gisc  (lie  output  spectra.  I  lie  program  can  be  quite  easily  extended  to  include  other  gust 
soniponenls:  cross-correlations  would  require  the  use  ol'  the  more  general  l  quation  (4)  for  the 
output  spectra. 


5.  COM  I  l  SIGNS 

Ibis  Note  has  described  a  I  ortrun  program  developed  for  aircraft  rigid-body  gust  response 
calculations.  I  he  theoretical  basis  and  Hs  limitations  base  been  summarised,  ss  it h  a  more  detailed 
description  ot  the  aircraft  and  gust  models  and  of  the  program  organisation. 

Decoupled  small  disturbance  equations  ol'  motion  base  been  used  ss uh  prosision  for 
additional  equations  for  any  (light  control  systems.  Most  Configuration.  I  light  Condition  and 
\crodsnumie  Data  are  read  in  as  inputs  and  can  readily  be  saried.  I  lie  gust  spectra  described 
are  appropriate  U>  isotropic  homogeneous  turbulence  and  cross-correlations  base  not  been 
included  in  the  examples  gisen. 

Ihe  main  program  has  been  organized  in  a  flexible  form  to  allow  easy  alteration  either  to 
the  mathematical  model  or  the  gust  spectra.  Thus,  for  example,  changes  in  flight  control  ss  stems 
oi  addition  of  extra  output  sariables  can  be  readils  accommodated.  Similarly  it  is  relatively 
c.iss  to  modify  the  program  to  cater  for  cross-correlations  in  the  gust  inputs.  The  calculation  of 
the  gust  response  transfer  functions  is  independent  and  separate  front  the  esaluation  of  the 
output  spectra  so  that  it  is  possible  to  use  the  program  simply  to  calculate  the  transfer  functions 
of  ans  linear  system. 

Ihe  examples  gisen  demonstrate  the  usefulness  of  the  program  in  assessing  the  effects  of 
variations  in  configuration,  aerodsnamics  or  control  systems  on  the  gust  response  of  a  gisen 
xehielc. 
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APPENDIX  1 


Aircraft  Dynamic  Model 


(i)  Basic  Model 

The  Longitudinal  A(.v)  matrix  is 
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1  he  lateral  At  v )  matrix  is 
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The  aerodynamic  transfer  functions.  6iv.  6'nv.  etc.,  arc  usually  replaced  by  a  "quasi- 
steady"  aerodynamic  deri\ati\e  representation,  i.e. 

(<t\  ~  (  ii 
(inv  =  Co  | 


but 


('  H/>  ~  f  ti,| 

etc. 


(/flu  -  (  Mr..  S(  in* 

(ii)  Model  with  Pitch  and  )'u\i  Pamper 

The  yaw  damper  equations  are  (Fqn  ( 1 1 )): 

•  ( t  ,  •  r:i).v  •  l  >«S(.v)  -  <S|l,v) 

( T I r:|.v-  •  (ti  r:I).v  ■  1  )»S , ( \ )  =  Ksq(s) 

The  state  vector  becomes  (Fqn  (b)): 

x  =  (A  F(v).  a(.v),  </(  v).  Af'(.y).  Sd.v),  .i.v)]r 

with  a  similar  extension  for  the  lateral  case  (F.qn  (7)). 


(A3) 


(A4) 


(Al) 


The  longitudinal  matrix  is  augmented  to 
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with  a  similar  extension  in  the  lateral  ease  except  that  (7 i»rf.  G'i„,  and  G’n«  are  replaced  hy 
Gjt.  G'u  and  G'nA  respectively. 


APPENDIX  2 


Input  to  Gust  Response  Program  [Mirage] 


(A)  Configuration  Data 

1 .  Mass,  m 

2.  Inertias 

(i)  Roll.  Is 

(ii)  Pitch,  /„ 

(iii)  Yaw,  I; 

(is )  Cross, 

3.  Reference  Area.  S 

4.  Reference  Chord,  f 

5.  Reference  Span,  A 

6.  Centre  Line  Chord 

7.  C.G.  Position  as  of  Centre  Line  Chord 
X.  Damper  Gain  (Lqn  (10)),  K 

(B)  f  light  Conditions 

9.  Airspeed,  V 

10.  Altitude,  h,  or  Density.  /> 

1 1.  Turbulence  Length  Scale.  L 

(Cl  Aerodynamics  Longitudinal 

12.  Co,,.  Trim  Drag 

13.  Cn, ,  Drag  due  to  speed 

14.  Cr>„.  Drag  due  to  incidence 

15.  Ci., .  Lift  due  to  speed 

16.  Ci.,.  Lift  due  to  incidence 

17.  Ci.tf.  Lift  due  to  pitch  rate 

IX.  C my.  Pitching  Moment  due  to  speed 

19.  Cm„.  Pitching  Moment  due  to  incidence 

20.  Cm;,,  Pitching  Moment  due  to  rate  of  change  of  incidence 

21.  Cmv.  Pitching  Moment  due  to  pitch  rate 

22.  Ci,s,  Lift  due  to  pitch  damper 

23.  Cms,  Pitching  Moment  due  to  pitch  damper 


(D)  Aerodynamics  -  Lateral 


24.  (\d.  Sideforce  due  to  sideslip 

25.  (  >(i.  Sideforce  due  to  roll  rate 

26.  Sideforce  due  to  yaw  rate 

27.  ( Rolling  Moment  due  to  sideslip 
2S.  Cir.  Rolling  Moment  due  to  roll  rate 
24.  C,,.  Rolling  Moment  due  to  yaw  rate 
50.  (  Rawing  Moment  due  to  slideslip 

31.  C  „;l.  Yawing  Moment  due  to  roll  rate 

32.  Ci,,.  Yawing  Moment  due  to  yaw  rate 

33.  (\ ,.  Sideforce  due  to  yaw  damper 

34.  (.  iv  Rolling  Moment  due  to  yaw  damper 

35.  c  5  awing  Moment  due  to  yaw  damper 


FIG.  1  BASIC  STRUCTURE  OF  'GUSTR' 
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FIG.  2  PITCH  RATE  RESPONSE  TO  VERTICAL  GUST,  aa  -  RPV 
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FIG.  3  PITCH  ATTITUDE  RESPONSE  TO  VERTICAL  GUST.  oB  -  RPV 


Log io  <“,» 


FIG.  4  ROLL  RATE  RESPONSE  TO  LATERAL  GUST,  AND  GUST  GRADIENT, 
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FIG.  5  ROLL  RATE  RESPONSE  TO  LATERAL  GUST,  0g  -  MIRAGE 
(YAW  DAMPER  OFF) 
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FIG.  6  YAW  RATE  RESPONSE  TO  LATERAL  GUST,  0g  -  MIRAGE 
(YAW  DAMPER  OFF) 
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FIG.  7  ROLL  RATE  RESPONSE  TO  LATERAL  GUST,  (3g  -  MIRAGE 

(Cn  =  -0.15) 


-3.0 


L°9'io 


frequency 

response 


Power  spectrum 
x  106 


Log10(U,) 


FIG.  8  YAW  RATE  RESPONSE  TO  LATERAL  GUST,  0  -  MIRAGE 

ICe  =  -0.1b) 
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